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Monte Carlo Study of the Separation of Energy Scales in Quantum 
Spin 1/2 Chains with Bond Disorder 



00 



C/3 



cd 



o 
o 



> 

in 
(N 
(N 

O 
On 
0^ 
-(— • 



O 

o 



X 



Beat AmmonIj and Manfred Sigrist 



ir*! 



SCSC and Theoretische Physik, Eidgenossische Technische Hochschule, CH-8092 Zurich, Switzerland 
^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502 

(Received ) 

One-dimensional Heisenberg spin 1/2 chains with random ferro- and antiferromagnetic bonds 
are realized in systems such as SraCuPti-xIrajOe- We have investigated numerically the thermo- 
dynamic properties of a generic random bond model and of a realistic model of SraCuPti-xIriOe 
by the quantum Monte Carlo loop algorithm. For the first time we demonstrate the separation 
into three different temperature regimes for the original Hamiltonian based on an exact treat- 
ment, especially we show that the intermediate temperature regime is well-defined and observable 
in both the specific heat and the magnetic susceptibility. The crossover between the regimes is 
indicated by peaks in the specific heat. The uniform magnetic susceptibility shows Curie-like be- 
havior in the high-, intermediate- and low-temperature regime, with different values of the Curie 
constant in each regime. We show that these regimes are overlapping in the realistic model and 
give numerical data for the analysis of experimental tests. 

KEYWORDS: quasi-one-dimensional spin system, random spin chain, quantum IVIonte Carlo, SrsCuPti—jjIrj^Oe 



§1. Introduction 

One-dimensional quantum spin cliains are typical ex- 
amples of many-body systems with a very rich variety 
of physical properties. Over many decades they have at- 
tracted much interest in theory and have motivated the 
development of various calculation schemes, both ana- 
lytical and numerical. The number of real compounds 
containing quasi-one-dimensional spin systems is grow- 
ing. AJjew examples are organic systems like NENP and 
NINO0& or inorganic compounds such as SrsMPtOe (M 
— Ni, Cu and Zn).EP Some of them are investigated as 
possible realizations of Haldane gap systems. Recently 
also ladder systems such as (Ca)2V03 and Sr2Cu03 were 
investigated as an example of a resonating valence bond 
system. All these systems are supposed to be regular 
chains or ladders. However, in practice disorder must 
occur in most of these systems. Already little disorder 
in the composition can have considerable influence on 
the low-energy properties. A peculiar example of a dis- 
ordered spin chain was recently discovered, by Nguyen 
and zur Loye: the alloy Sr3CuPti_a;Ira;06£l-' The pure 
compound SraCuPtOg forms an antiferromagnetic (AF) 
spin chain where a spin 1/2 is provided by each Cu-ion. 
In this compound the Cu-ions are alternating with the 
spinless Pt-ions along chains. If all Pt is replaced by 
Ir which carries a spin 1/2, then the system is a fer- 
romagnetic (FM) spin chain. Therefore the random al- 
loy Sr3CuPti_a:Ir3:06 is a spin system with two types 
of bonds, FM (Jp < 0) and AF (Ja > 0), which are 
randomly distributed. There is a correlation among the 
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bonds in the sense that FM bonds always occur in se- 
quences of even numbers, since each Ir-ion makes such 
bonds with its two neighboring Cu-ions. A generic model 
to study the properties of such a system was considered 
by Furusaki and coworkers^ who neglected the correla- 
tion among the bonds and used J = |Jp| = |Ja|- The 
quantum spin system in their analysis has the following 
nearest neighbor Heisenberg Hamiltonian 



'H — 2_^ Ji^i ■ Si+l, 



(1.1) 



with a bond probability distribution which is indepen- 
dent of i 

P{.h) = p5{J, + Jf) + (1 - p)5{.h - Ja), (1-2) 

where < p < 1, and 5 is the Kronecker delta function. 
For p = Q and 1 we have a purely AF or FM spin chain, 
respectively. 

The ajialysis of this model by high-temperature ex- 
pansioncr and W a real space renormalization group 
(RSRG) schemeB' suggests that three different temper- 
ature regimes are present in this system. In the high- 
temperature regime the spins essentially behave as in- 
dependent degrees of freedom. If the temperature is 
lowered down to /cgT ^ J the spins start to correlate. 
First they align within the segments of purely FM or AF 
bonds. The coupling among the segments is weak. Thus, 
the spins in each segment create an effective single spin 
degree of freedom which is rather large in the case of FM 
bonds and 5 = or 1/2 for AF bonds. At intermediate 
temperatures these effective spins behave as independent 
due to the thermal fluctuations. However, they begin to 
correlate at lower temperatures. The intermediate and 
low-temperature regimes are described by a random spin 
model with a basically continuous random bond distri- 
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bution as well as random spin sizes (effective spin sizes). 

The limit of T ^ was analyzed by Westerberg and 
coworkeiis by means of a RSRG scheme for this type of 
modelfl'EP Here we will briefly review their main results. 
As the temperature T is lowered, a growing number of 
spins is correlated in clusters consisting of n spins on the 
average, which form a large spin S. The spin size S scales 
with S ex n}'"^ and the average number of spins in a clus- 
ter scales with the temperature as ri oc T~^". ThcLscaling 
exponent a « 0.22 was determined numerically .BlP Be- 
cause the large effective spins S of the clusters behave as 
essentially independent under thermal fluctuations, one 
expects a Curie-like susceptibility x ^ C/T with a T 
dependent Curie constant C. 

From this result we conclude that, in principle, this 
system exhibits three different temperature regimes, each 
with its own Curie-like susceptibility and Curie constant. 
In the high-temperature regime {k^T 3> J) the Curie- 
behavior originates from the S = 1/2 spins leading to 

X = f^: where /ig = 27^ ^^ ^^^ Bohr-magneton and 
the Curie constant is C = S{S + l)/3 == 1/4. In the in- 
termediate temperature regime {ksT ^ J) the effective 
spins of the segments give a p-dependent Curie constant 
Cgff. Finally a crossover to a scaling regime with a new 
Curie behavior occurs at very low-temperatures. Obvi- 
ously the three Curie constants should be a decreasing 
sequence as the number of available degrees of freedom 
is decreasing with lowering temperature. 

The separation into three regimes is expected to be 
visible in the specific heat too. Peak- or shoulder-like 
structures indicate the boundaries between the regimes 
as they are sign of correlation of degrees of freedom. One 
boundary occurs near k^T ^ J and the second at some 
lower temperature. In the scaling regime at very low 
temperatures, the assumption of independent large spins 
S of correlated clusters leads to cy/T oc T^°'-^\\tlT\Q> 

The crossover between the high- and intermediate 
temperature regime has been examined by the high- 
temperature expansion. For low temperatures, however, 
only the limiting behavior for T ^ has been inves- 
tigated for an effective Hamiltonian-with a broad ran- 
dom distribution of couplings JiP'tP which is in con- 
trast to the discrete distribution of the initial spin cou- 
plings in Sr3CuPti_a:Ira:06. From these results we know 
that the effective spin scaling regime starts at very low 
temperatures, which are hard to observe by experiment. 
The actual crossover from the intermediate to the low- 
temperature regime can be investigated more easily by 
experiment, but has not yet been analyzed so far. In this 
paper we would like to close this gap by demonstrating 
that the intermediate temperature regime is well-defined 
and observable in both the specific heat and the suscepti- 
bility by an exact treatment of the original Hamiltonian. 

The correlations among the FM bonds in 
Sr3CuPti_j;Irj;06 are another important aspect which 
has not yet been examined. We will investigate the 
implications of these correlations in a realistic model 
where the FM bonds always occur pairwise and 
study its differences to the uncorrelated model in the 
magnetic susceptibility and the specific heat in an 



experimentally accessible temperature range. These 
data will provide a sensitive test for the randomness of 
the distribution of the Pt-ions in SrsCuPti^ajIrajOg. 

Our investigations are based on QuantunxMonte Carlo 
(QMC) simulations by the loop algorithm,© which is an 
excellent method to simulate accurately the thermody- 
namics of large spin systems. 

§2. Numerical Methods 

The QMC loop algorithmcl^ is a finite temperatuxe 
method based on the Trotter-Suzuki decompositiont3^ 
of the partition function Z and is ideally suited for the 
calculation of thermodynamic properties of unfrustrated 
spin systems. It allows the direct and exact compu- 
tation of various thermodynamic observables, e.g. the 
uniform magnetic susceptibility x or the internal energy 
u, without introducing any approximations. In contrast 
to the classical Metropolis world-line algorithmjEJ' the 
loop-algorithm does not suffer from prohibitively long 
auto-correlation times, j— Additionally, the introduction 
of improved estimatoral3) gives a further reduction of 
the variance. This enables us to investigate much bigger 
problems than with previous QMC methods. 

For the calculation of the random bond models, we 
have considered up to 400 samples of different random 
bond configurations, each consisting of a chain of i = 
400 sites with periodic boundary conditions. In Fig. 1 
we show the distribution of FM (AF) segments of length 
Is in our configuration samples. We have performed 
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Fig. 1. Statistic of the total distribution of the segments with Is 
FM (squares) and AF (circles) bonds in the configurations we 
have used for the simulation of the unconstrained model. The 
dashed curve shows the ideal distribution. 



calculations in the temperature range from T/J — 1/60 
to T/J = 10. Unless otherwise mentioned, we have used 
Trotter numbers between 20 and 120 and extrapolated to 
Trotter time-steps At J = 0. We made 10^ multi-cluster 
updates for thermalization, followed by 2 x 10^ up to 2 x 
10® multi-cluster updates for measurements, depending 
on the temperature T and the Trotter-number. 
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The value of an observable O is estimated by averag- 
ing over M successive measurements Oj coming from the 
QMC simulation: (O) « O = M"! YljLi Oj, and the er- 
ror is determined by the variance of the measurements. 
In this way we have measured directly the internal energy 
u, the total magnetization M, and the uniform suscep- 
tibility X for each configuration. Due to the efficiency 
of our QMC algorithm, we have been able to measure 
the internal energy u accurately enough to calculate the 
specific heat numerically as the first derivative of the in- 
ternal energy u with respect to the temperature T . After 
the calculation of all observables for each individual con- 
figuration, the average over all configurations has been 
taken. In order to get a correct estimate of the errors, we 
have taken into account the variance of the observable for 
the set of configurations and the error of the observable 
for a single configuration. In general, the error-bars for 
the observables of a single configuration are much smaller 
than the variance of the observables for the set of config- 
urations, and therefore we can neglect them completely. 

§3. Results 

We have studied two different models. First, we have 
investigated a generic random bond model with Ti. = 
'Y^^ JiSi ■ Si+i and the bond probability distribution 

P{J,) = p6{.h + Jf) + (1 - P)&{'h - J a). (3.1) 

where < p < 1 and |JfHI>^aI for the case of p = 0.5, 
i.e., the same probability for finding a FM as an AF bond. 
This model has beerLsludied in a number of papers using 
different methods Jj'BEP In the rest of this paper we will 
call this model the "unconstrained" model. 

In the real Sr3CuPti_a;Ira;06 alloys, there are some ad- 
ditional restrictions. As we have previously mentioned in 
the introduction, by replacing the spinless Pt-ions with 
Ir-ions, two FM bonds are created between the spin 1/2 
carrying Ir-ions and its neighboring Cu-ions. In addition 
to this pairwise correlation of FM bonds, the FM cou- 
plings are stronger than the AF couplings, an estimate 
gives I JfI = 4| Ja|. By generating a random sequence of 
Pt- an Ir-ions we have generated a configurations equiva- 
lent to and Ir substitution of a; = 0.5, which corresponds 
to the probability of finding a FM bond with p = 2/3. 
We will refer to these systems as the "constrained" model 
in the following text of this paper. 

§4. Susceptibility 

Let us consider the unconstrained model first. In the 
high-temperature regime /cgT ^ J, we have indepen- 
dent single spins as in the uniform AF or FM chain, and 
a Curie-law for the susceptibility Xc\{T) = M^/(4fcB2^) for 
T —>■ oo. As the temperature T is lowered, the individual 
spins start to correlate. In the AF segments, they form 
the collectively lowest possible total spin {S = for an 
odd number of bonds Is oi S = 1/2 for an even number 
Is), while the FM segments form the collectively largest 
spin (by aligning all spins parallel). Due to the mis- 
fit of the discrete spectra in the AF/FM spin segments 
(the finite size gap SE-ps is 5Epg oc J/ls for AF and 
^E'fs °^ J/^l foi' F^'l segments), the excitations remain 



localized in the segments and the interactions among dif- 
ferent segments remain very weak. In this intermediate 
temperature regime the segments essentially behave as 
uncoupled effective spins 5'eff. These uncoupled effective 
spins can be seen in a second Curie-law of the suscepti- 
bility X on an intermediate energy-scale (ex SEpa,) with 
an effective Curie-constant Cgg depending on the average 
size of the effective spins 



X 



knT 



with Ceff = i g™it£) '^ ^he^e (%M/Af) denotes 
the effective spin of the FM/AF segment, and (fiFM/Ap) 
is the average length of a FM/AF segment. The spins 
on the boundary of a segment of pairs of FM spins and 
a segment of AF spins always have to be counted to 
the AF segment and not to the FM segment, i.e., the 
total spin of a FM segment with a length of Is bonds is 
'5'tot — S{ls — l).tP For the unconstrained model with 
p — 0.5 and S = 1/2 one obtains C^q = 1/8 for a chain 
of infinite length. 




k,T/J, 



Fig. 2. Figure of the uniform magnetic susceptibility x times 
temperature T. The symbols are the results of the QMC simula- 
tions for the unconstrained model (squares) and the constrained 
model of SrsCuPti-^ilri^Oe (diamonds). For comparison we 
also show the results for the uniform FM and AF Heisenberg 
chain. The Curie-law behavior can be seen by horizontal lines, 
the dashed lines show the results of the statistical analysis, sq 
are the limits of x^ for T ^ 0. 



In Fig. 2 we show xJ"^ *s a function of the temper- 
ature, such that the Curie-law behavior of the effective 
spins can be seen as a plateau. From our QMC simula- 
tions we get Cgff = 0.138 ± 0.003, whereas the value of 
the effective constant Ceff obtained by statistical anal- 
ysis of the bond distribution is C^g = 0.1252, and it is 
shown by a dashed line in the graph. Previous results 
by the transfer matrix method apd high-temperature ex- 
pansions obtained C^g = 0.13.EP However, in the high- 
temperature expansion, a possible error stems from the 
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extrapolation of the Pade approximants. There is also 
some range for modification by considering to which seg- 
ment a spin on the border between FM and AF segments 
has to be counted in the statistical analysis: In general, 
the spins at the border between these segments tend to 
be included in the AF segments, but if a long AF segment 
is next to a short FM segment, some of these boundary 
spins rather couple to the FM segment, resulting in an 
increase of C^s- A detailed analysis of the bond distri- 
bution used in our QMC calculations shows that if we 
count the spins on the border of AF segments consisting 
of six and more bonds next to FM segments consisting of 
three and less bonds to the FM segment, we rather get 
Cgff = 0.131 and similarly for AF segments consisting 
of five and more bonds next to FM segments consisting 
of three and less bonds we find Cgg — 0.138 instead, in 
agreement with our QMC results. 

The temperature range of the regime of uncoupled 
effective spins is determined by the extension of this 
Curie-law plateau, it starts at k-gT « 0.16 J and ends 
at kj^T « 0.06J. The region of this regime compares 
favorably to the simple estimate given by the finite size 
gap of the spin segments SEpg . 

At very low temperatures, interactions among the ef- 
fective spins become relevant and the effective spins start 
to correlate. These interactions are again described by a 
Heisenberg model with H^g = ^^ •^eff^eff ' ^eff^' "^^^re 
the sum i is now taken over the spin segments. The 
couplings J'rr are random in both sign and magnitude, 
and the size of the effective spins S^^ff is random as well. 
Following ref. 6, the distribution of the effective spins 
S'V and their interactions J*rr is very broad. Hence 
two neighboring spins with the largest energy gap to 
their first excited state will be locked into their ground- 
state, consisting of their maximal (Jgff < 0) or min- 
imal (JeS > 0) spin state and form a new effective 
spin Sf,g. As more and more spins freeze out, the ef- 
fective spin size S^g scales with the number of origi- 
nal spins in the cluster as S'gg ~ n^^"^ (random walk). 
The average number of original spins in a cluster n de- 
pends on the temperature T and scales with a power law 
n ~ T-2", with a = 0.22 ± O.OlJj) consistent with the 
resuh a = 0.21 ± 0.01 of ref. 8 

In Fig. 3 we have plotted JxT as a function of T°'/J 
for different values of a. The low temperature behavior 
has been calculated in ref. 8: 



1 



(4.1) 



where sq is the groundstate expectation value of the total 
spin per site 
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M 



^i—l 



S,Si 



(4.2) 



where Si^i = —SiSgn{Ji) and Si = I. The beginning of 
the low-temperature scaling regime is estimated to start 
at T < 0.05 Jo for the effective Hamiltonian,® where Jq 
is the maximum of a broad random distribution of cou- 
plings. This distribution of couplings corresponds to the 
effective couplings Jefi of the intermediate temperature 
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Fig. 3. Low temperature scaling behavior of x^ for the un- 
constrained model (a) and the constrained model (b). The low- 
temperature scaling regime is characterized by x'^ = sq+T" /'y+ 
0{T^°'). The open symbols denote different values of the scaling 
exponent a, the full symbols denote the T ^ limit 



regime. These couplings are about an order of magni- 
tude smaller than the original couplings Jj^ In Fig. 3, 
we can see that a linear extrapolation of the data points 
with the scaling exponent a ~ 0.22 does not yet yield the 
correct zero-temperature value sq and we conclude that 
the scaling regime has not been reached for T = 0.02J. 
Unfortunately, critical slowing down in the QMC sim- 
ulation at low temperatures prevents us from reaching 
temperatures as low as T < 0.005 J. 

Next we concentrate on the susceptibility of the con- 
strained model of Sr3CuPti_a:Ira;06. In this case we have 
to take into account the different strength of the FM and 
AF interactions. Starting from the same high tempera- 
ture limit x/L — jJ/T as in the generic case, the FM 
bonds start to correlate by aligning parallel at higher 
temperatures T than the AF bonds, due to the stronger 
coupling |JfI = 4|Ja| = 4J. The formation of these 
large spins leads to an increase of JxT first. Then the 
AF spins start to correlate too, but instead of a plateau 
as in the case of the unconstrained model we obtain a 
broad peak near kj^T = J. 

If we calculate the average effective spins {S-pyi/j^-p) 
as in the previous case, we obtain an estimate of the 
effective Curie constant Cgff = 0.33 ± 0.017. This 
value is a lower bound of the Curie constant, where 
the spins within the segments are completely correlated. 
The effect of the constraint among the FM bonds in 
Sr3CuPti_2,Ira;06 is very important, if we neglect that 
we would obtain Cgg « 0.2731 instead. An upper bound 
of C^^ = 0.40 ± 0.02 can be obtained by assuming 
that only the spins in the FM segments correlate and the 
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AF segments remain uncoupled due to their weaker cou- 
plings. From our QMC data in Fig. 3 we get C^s ~ 0.36. 
The deviation from C^g of the spin segment statistic is 
larger than for the unconstrained model, since the dif- 
ference in the coupling strength favors quantum effects 
on the boundaries of large FM and small AF segments 
and makes the simple statistical estimate more subtle. 
The effects causing an increase of of Cgfj are twofold in 
this case: the first is the same as for the unconstrained 
model, spins on the border of large AF segments next 
to a short FM segments tend to couple to the FM seg- 
ment. A detailed analysis of the spin distribution gives 
C'eff — 0-34 if we count the spins on the border of AF 
segments consisting of five and more bonds next to FM 
segments consisting of four and less bonds to the FM 
segment, and C^fj ~ 0.345 for AF segments consisting 
of four and more bonds next to a segments consisting 
of three and less bonds. The second reason for a larger 
value of Ccfj is that due to the relatively weaker AF cou- 
plings than the FM couplings, some of the AF segments 
are not yet completely correlated. Thus especially the 
long AF segments still consist of some original spin-1/2 
degrees of freedom, giving rise to a further increase of 
Ceff as we have demonstrated for the case of completely 
uncorrelated AF segments before. 

Additionally we have also determined Cgfj by a QMC 
simulation for the same distribution as above but with 
different coupling strength \Jp\ = 2|JaI = 2 J. Here we 
can see again how the susceptibility x'l'J increases from 
the high-temperature limit x'l'J = 1/4; but in contrast 
to the previous case it never reaches the lower bound of 
Cgff = 0.33 ± 0.017 of the regime of uncoupled effective 
spins. In this case, some spin segments already start to 
correlate and form a collective spin before the regime of 
independent effective spins is established. 

In order to investigate the low-temperature scaling 
behavior of eq. (4T) for the cons trained model, we 
have to calculate sq from eq. (42). Let us consider 
a chain of L Cu- ions. There are Lx spin 1/2 car- 
rying Ir-ions in the chain and the ground state ex- 
pectation value of the total spin is: ((X^spins '^'i)^) — 
i^ + x)L{Sf)+J2spms(^i^j)' where the sum is taken over 
all Cu- and Ir-ions. For the evaluation of X]spins('^j'^j) 
we divide the sum into four terms using the follow- 
ing notation: we enumerate all Cu-ions with the in- 
dex i, while we use the index i for the Ir-ion between 
the Cu-ions i and i + 1 and write J^sp'msi^i^j) ~ 



'^i^jiSiS-Aj. Now we can calculate the expectation 



value (SiSj) by taking into account that an Ir-ion exists 
only with the probability x and that each missing Ir-ion 
changes the sign of the next spin: 



{S^S, 



i+nl 



\Y,x"^{i^xr-^{-i) 



n — rn/ i \7i—m 



-{2x-ir 



(^*^.rJ = jE^'"(i-^)""™(-i)' 
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'-{2x-iy 



n-l 



(4.3b) 



1 " " 
{SiS,+n) = -^x™(l-a:)"-i-'"(-l)"-i- 



m=0 



n- 1 
m 



-{2x-ir'' 



(4.3c) 



'^^l^Tr) = f ^2:"(l-x)"-i-™(-l)"-i- 



m 



-{2x-ir-\ 



(4.3d) 



where ( '" ) — n\ / {in\{n — m)\) . From this we obtain 



^5,1 )^L{l+p){S^) 

Y,{'^x-iT + Y,x{2x-ir 

.71—1 n— 

+x [ X^(2x - l)"-i -f ^ x{2x - l)"-i j 
= L, (4.4) 



^ spins 

2L 

T 



groundstate 
0.222222.... 



which gives in the 

■•^O = ((Espins^0')/(3fcB(l + ^)i) = 2/9 

In Fig. 3 we can see that also for this model the value 
we get for sq is too small if we extrapolate our data lin- 
early with the correct scaling exponent a ~ 0.22. Thus 
we conclude that the low-temperature scaling regime 
starts below T/Jj^ < 1/100. 

§5. Specific Heat 

Again, let us concentrate on the unconstrained model 
first. In the specific heat data, the crossover between 
different temperature regimes should be visible by peak- 
like structures, where the spins start to correlate. In the 
uniform FM or AF Heisenberg systems, one broad peak 
below k^T/J ~ 1 appears, where the individual spins 
start to correlate. Two peaks are expected in the ran- 
dom system. One where the original 5 = 1/2 spins start 
to correlate (as for the uniform systems), and a second 
one where the segments of the effective spins start to cor- 
relate (at the end of the Curie-plateau k-Q,T/ J ^ 1/20). 
In Fig. 4 we show the specific heat per spin cy of the ran- 
dom bond systems and of the uniform AF and FM chain 
for comparison. The large and broad peak is the sig- 
nature of the correlation of the individual spins. In the 
inset, one can clearly see the second peak at low tem- 
peratures of the correlation of the effective spins for the 
unconstrained model. This confirms the clear separation 
of the different regimes and could not be investigated 
(4.3a) with the previous methods, e.g. the high-temperature 
expansion. 

In the scaling regime at very low temperatures T, 
the assumption of independent large spins S of cor- 
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Fig. 4. The specific heat per spin cy of the unconstrained ran- 
dom bond model and the SraCuPti-^^IriiOe-model. For com- 
parison, we have also shown the results of the uniform FM and 
AF Heisenberg chain. The symbols are the results of the QMC 
simulation for the generic model (squares) and the constrained 
model of SrsCuPti-a^Ira^Oe (diamonds), dashed and dotted lines 
(HTE) are the high-temperature expansion results obtoined by 
the extrapolation of the two-point Fade approximants.Q In the 
inset, one can see the peak in the specific heat of the uncon- 
strained model and a cusp for the constrained model which in- 
dicate the onset of correlations among the effective spins. 



related clustersB leads to the entropy per site a oc 
ln(25'(T) + l)/n{T) oc T-2"|lnr|. From this one finds 
cy/T (X T^^"^^|lnT|. But as for the magnetic suscep- 
tibility, we can not reach such low temperatures as to 
determine the scaling exponent a from our data. By the 
relation /„ dTcy/T, we can see from the divergence of 
cy/T as T — > that a large fraction of the entropy is 
at very low temperatures, due to the broad spectrum of 
energy-scales of the effective spins. In Fig. 5 we show 
cy/T, the area below the curve corresponds to the en- 
tropy S. If we consider the data in Fig. 5 in more detail, 
we can see that the value of cy/T for the generic model 
is always smaller than cy/T of the uniform AF Heisen- 
berg chain. Only at very low temperatures, we can see a 
sharp crossover to the divergent low-temperature behav- 
ior. This crossover is due to the very different energy- 
scales of the effective spins and can not be seen by the 
high-temperature expansion. 

More important for comparisons with experiments 
is the specific heat of the constrained model of 
SrsCuPti-ajIrajOe. Here, the spins start to correlate at 
different temperatures T/J due to the different coupling 
strengths. Hence there is no clear transition and the peak 
at T/J « 1 is much weaker than for the unconstrained 
model. Also the peak at the transition to the scaling 
regime is much weaker, but we can interpreted the cusp- 
like structure near T « 0.05 J as the corresponding peak 
in the specific heat. In Fig. 5, where we show cy/T, one 
finds that there is no peak in cy/T before the crossover 



Fig. 5. The specific heat per spin cy divided by the temperature 
T for the uniform and random bond systems. The area below 
the curves corresponds to the entropy S. The symbols are the 
QMC results, the dot-dashed line (HTE) is the high-temperature 
expansion result obtained by the extrapolation of the two-point 
Fade approximants.U^ One can see that the QMC and high- 
temperature expansion results are in perfect agreement down to 
fcgT Ki 0.2J, then the QMC results show the correct crossover 
to a divergent behavior. 



to the divergent low-temperature behavior as for the un- 
constrained model, but cy/T is rather continuously in- 
creasing for r -^ 0. 

§6. Conclusions 

We have investigated numerically the thermodynamic 
properties of one dimensional spin chains with random 
FM or AF couplings. We have considered an uncon- 
strained model with equal magnitude of the AF and 
FM couplings and a realistic model of Sr3CuPti_a:Ir2;06 
(constrained model). The thermodynamic properties of 
the random bond systems are determined by three differ- 
ent energy scales: at high temperatures, the energy scale 
is given by the original spin-spin couplings. The suscep- 
tibility obeys a Curie law of free spins at very high tem- 
peratures (T ^ k^T). By lowering the temperature, the 
spins start to correlate within the AF and FM segments 
and form effective spins. At intermediate temperatures, 
these effective spins interact very weakly, essentially as 
free spins. The crossover of the individual spins to effec- 
tive spin segments can be seen by a peak in the specific 
heat and a Curie-law of the magnetic susceptibility of 
the effective spins. The energy scale of this intermediate 
temperature regime is given by the interactions of these 
effective spins. At still lower temperatures, the inter- 
actions among the effective spins become relevant and 
the spin segments gradually freeze out into clusters of 
correlated spin segments. Also this crossover to the low- 
temperature scaling regime is visible by a peak in the 
specific heat, although this peak is much smaller than 
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for the previous transition. 

Our results represent the first exact treatment of 
the original Hamiltonian showing a clear separation 
into three different temperature regimes for the uncon- 
strained model. This separation can be seen by two 
marked peaks in the specific heat and two Curie laws 
in the uniform magnetic susceptibility at high and in- 
termediate temperatures. In order to allow the compari- 
son with experimental results, we investigate for the first 
time a realistic model of SrsCuPti^jjIrajOg with a con- 
strained distribution of the FM bonds exactly. 

We find that the three energy levels are not so clearly 
separated in this case. Because of the different magni- 
tude of the FM and AF couplings, the energy scales are 
overlapping and the crossover between the three regimes 
is continuous. Therefore the peak originating from the 
correlation of the original spins is very broad and the 
onset of correlations among the effective spins can only 
be seen by a small cusp in the specific heat. Instead of 
a clear plateau in xT ^ the effective spins form a broad 
peak. 

Overlapping energy-scales can be found in many ran- 
dom bond systems, indeed, the formation of a clearly dis- 
tinguishable intermediate temperature regime depends 
subtly on the distribution and strength of the FM/AF 
bonds. However, the general picture still applies also 
for systems with overlapping energy scales and we have 
given upper and lower bounds for the Curie-constant of 
the susceptibility in the intermediate temperature regime 
of the constrained model. The calculation of these Curie- 
constants can be used as a universal tool for the anal- 
ysis of the bond-distribution, e.g. much larger Curie- 
constants will result in experiments if the distribution of 
Pt- or Ir-ions is not really random because of large FM 
clusters which can enhance the Curie-constant drasti- 
cally at low temperature. The calculation of these Curie- 
constants can also be appliedJii_.the statistical cluster- 
analysis of bond-distributions .BI13-' 

We have given upper bounds for the beginning of 
the low-temperature scaling regime, T = 0.02J for 
the generic model and T — O.OUa for the realistic 
model. But despite the very efficient QMC algorithm, we 
have not been able to reach the low-temperature scaling 
regime with the original Hamiltonian for both models. 
For the original Hamiltonian, the scaling regime starts 
at extremely low temperatures, such that a completely 
different approach will be needed to reach this regime 
with the original Hamiltonian. However, our numerical 
results confirm that the intermediate temperature regime 
is well described by an effective Hamiltonian for the ef- 
fective spins, providing a firm basis for the investigations 
of the scaling regime based on the effective Hamiltonian. 

We would like to finish by mentioning that an analogue 
of the random bond spin chains exists £01 very dilutely 
randomly depleted Heisenberg ladders .O* Due to the 
spin gap in the parent material, the Curie law of the ef- 
fective spins can be observed experimentally more easily 
and the value of the Curie constant is in agreement with 
the theoretical valucEj* 
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